Steering quantum transitions between three crossing energy levels 
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We calculate the propagator and the transition probabilities for a coherently driven three-state 
quantum system. The energies of the three states change linearly in time, whereas the interactions 
between them are pulse-shaped. We derive a highly accurate analytic approximation by assuming 
independent pairwise Landau-Zener transitions occurring instantly at the relevant avoided cross- 
ings, and adiabatic evolution elsewhere. Quantum interferences are identified, which occur due to 
different possible evolution paths in Hilbert space between an initial and a final state. A detailed 
comparison with numerical results for Gaussian-shaped pulses demonstrates a remarkable accuracy 
of the analytic approximation. We use the analytic results to derive estimates for the half-width of 
the excitation profile, and for the parameters required for creation of a maximally coherent super- 
position of the three states. These results are of potential interest in ladder climbing in alkali atoms 
by chirped laser pulses, in quantum rotors, in transitions between Zeeman sublevels of a J = 1 level 
in a magnetic field, and in control of entanglement of a pair of spin-1/2 particles. The results for 
the three-state system can be generalized, without essential difficulties, to higher dimensions. 

PACS numbers: 32.80.Xx, 33.80.Be, 32.80.Rm, 33.80.Rv 
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I. INTRODUCTION 

Whenever the energies of two discrete quantum states 
cross when plotted against some parameter, e.g. time, 
the transition probability is traditionally estimated by 
the famous Landau-Zener (LZ) formula [1]. Although 
the LZ model involves the simplest nontrivial time de- 
pendence - linearly changing energies and a constant in- 
teraction of infinite duration, when applied to real phys- 
ical systems with more sophisticated time dependences 
the LZ model often provides more accurate results than 
expected. This feature (which has not been fully under- 
stood yet), and the extreme simplicity of the LZ transi- 
tion probability, have determined the vast popularity of 
the LZ model, despite the availability of more sophis- 
ticated exactly soluble level-crossing models, e.g. the 
Demkov-Kunike model [2] and its special case, the AUen- 
Eberly-Hioe model [SJ. 

Numerous extensions of the LZ model to multiple levels 
have been proposed. The exactly soluble multistate LZ 
models belong to two main types: single-crossing bow-tie 
models and multiple-crossings grid models. In the bow- 
tie models, where all energies cross at the same instant 
of time, analytic solutions have been found for three Q 
and N states [a, H 0] ' ^^'^ when one of the levels is split 
into two parallel levels @] . In the grid models, a set of Na 
parallel equidistant linear energies cross another set of N}, 
such energies (Dcmkov-Ostrovsky model) ^9i. I lOl [ill [l3| . 
For Nf, — I (or Na — 1) the Demkov-Ostrovsky model 
reduces to the Demkov-Osherov model [13, [ij]- The cases 
of one [1^ or two [la| degenerate levels have also been 
solved. In the most general case of linear energies of 
arbitrary slopes, the general solution is not known, but 
exact results for some survival probabilities have been 
derived [13, [H, [S [13] ■ 

A variety of physical systems provide examples of mul- 



tiple level crossings. Among them we mention ladder 
climbing of atomic and molecular states by chirped laser 
pulses I2l[,[22|, harpoon model for reactive scattering [23| . 
and optical shielding in cold atomic collisions [24!|. Ex- 
amples of bow-tic linkages occur, for instance, in a rf- 
pulse controlled Bose-Einstein condensate output coupler 
|25l [26| and in the coupling pattern of Rydberg sublevels 
in a magnetic field ^61]. A degenerate LZ model emerges 
when the transition between two atomic levels of angular 
momenta Ja and Jb = Jq or Ja ± 1 is driven by linearly 
chirped laser fields of arbitrary polarizations [1^ [l^| . 

A general feature of all soluble nondegenerate mul- 
tilevel crossing models is that each transition proba- 
bility Prn^n between states tpm and tpn is given by a 
very simple expression, as in the original LZ model, al- 
though the derivations are not trivial. In the grid mod- 
els, in particular, the exact probabilities Pm^n have the 
same form (products of LZ probabilities for transition or 
no-transition applied at the relevant crossings) as what 
would be obtained by naive multiplication of LZ proba- 
bilities while moving across the grid of crossings from ipm 
to ipn, without accounting for phases and interferences. 
Quite surprisingly, interferences between different paths 
to the same final state, a multitude of which exist in the 
grid models, are not visible in the final probabilities. 

In this paper we develop an analytic description of a 
three-state model wherein the three energies change lin- 
early in time, with distinct slopes, thus creating three 
separate level crossings. This system is particularly con- 
venient for it presents the opportunity to investigate 
quantum interference through different evolution paths 
to the same final state, and in the same time, it is suf- 
ficiently simple to allow for an (approximate) analytic 
treatment; for the latter we use sequential two-state LZ 
and adiabatic-foUowing propagators. This system is also 
of practical significance for it occurs in various physical 



situations, for instance, in transitions between magnetic 
sublevels of a J = 1 level [2^, in chirped-pulse ladder 
climbing of alkali atoms [23 , in rotational ladder climb- 
ing in molecules l28ll. and in entanglement of a pair of 
spin-1/2 particles [2^. The results provide analytic esti- 
mates of all nine transition probabilities in this system. 
We do establish quantum interferences and estimate the 
amplitude and the frequency of the ensuing oscillation 
fringes, as well as the conditions for their appearance. 
The analytic results also allow us to prescribe explicit 
recipes for quantum state engineering, for example, to 
create an equal, maximally coherent superposition of the 
three states. 

This paper is organized as follows. In Sec. |ll]we pro- 
vide the basic equations and definitions and define the 
problem. In Sec. Illll we derive the propagator, the tran- 
sition probabilities and the validity conditions. In Sec. 
IIVI we compare our analytical approximation to numeri- 
cal simulations. Then in Sec. |V]we demonstrate various 
applications of the analytics. In Sec. IVII we compare 
our model with the exactly soluble CarroU-Hioe bowtie 
model in the limit of vanishing static detuning. Finally, 
we discuss the conclusions in Sec. I VIII 



II. DEFINITION OF THE PROBLEM 

A. Description of the system 

We consider a three-state system driven coherently by 
a pulsed external field, with the rotating-wave approxi- 
mation (RWA) Hamiltonian (in units h — 1) 
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Ao + At ifii2(i) 
ir!i2(t) iriaaW 

ln23{t) Ao-At 



(1) 



The diagonal elements are the (diabatic) energies (in 
units h) of the three states, the second of which is taken 
as the zero reference point without loss of generality. Aq 
is a static detuning, and ±At are the linearly changing 
terms. To be specific, we shall use the language of laser- 
atom interactions, where the difference between each pair 
of diagonal elements is the detuning for the respective 
transition: the offset of the laser carrier frequency from 
the Bohr transition frequency. The pulse-shaped func- 
tions Qi2{t) and fl23{t) are the Rabi frequencies, which 
quantify the field-induced interactions between each pair 
of adjacent states, ■0i ^^ "02 and ^2 <-^ V^Sj respectively. 
Each of the Rabi frequencies is proportional to the re- 
spective transition dipole moment and the laser electric- 
field envelope. As evident from the zeroes in the corners 
of the Hamiltonian ([T]) we assume that the direct transi- 
tion tpi <-> ip^ is forbidden, as it occurs in free atoms when 
ipi <-^ "02 and 'ip2 ■*-> V'3 are electric-dipole transitions. 
The probability amplitudes of our system C(i) = 



c 
w 




Time (arb. units) 

FIG. 1: (Color online) Diabatic and adiabatic energies vs 
time for a Gaussian-shaped laser pulse. The labels denote 
the respective diabatic and adiabatic states. 



[Ci (i) , C2 (t) , C3 (i)] satisfy the Schrodinger equation 



iC{t) = H(i)C(t), 



(2) 



where the overdot denotes a time derivative. 

Without loss of generality, the couplings 1^12 (i) and 
^23(0 are assumed real and positive and, for the sake of 
simplicity, with the same time dependence. For the time 
being the detuning Aq and the slope A are assumed to 
be also positive. 



Ao > , ^ > 0; 



(3) 



we shall consider the cases of negative Aq and A later 
on. With the assumptions above, the crossing between 
the diabatic energies of states ipi and "02 occurs at time 
t_ = — T, where t — Aq/A, between ^2 and ip^ at time 
t+ = T, and the one between ^i and ips at time to = 0. 

Fig. [1] plots diabatic and adiabatic energies vs time 
for a Gaussian-shaped laser pulse. We use ipk and Xfc to 
denote diabatic and adiabatic states, respectively. 

The objective of this paper is to find analytical ex- 
pressions for the evolution matrix and for the transition 
probabilities between different diabatic states. 



B. Implementation 

The Hamiltonian ([1]) appears naturally in a number of 
specific problems of interest in time-dependent quantum 
dynamics of simple systems. 

The first example is ladder climbing of electronic en- 
ergy states in some alkali atoms, for instance, in rubid- 
ium [271 . A linearly chirped laser pulse couples simul- 
taneously both transitions 5s-5p and 5p-6s. If the car- 
rier frequency of the pulse is tuned on two-photon res- 
onance with the 5s-6s transition, then the intermediate 



state 5p remains off resonance, by a detuning A, which 
leads to the "triangle" Unkage pattern in Fig. [T] The cou- 
plings rii2(i) and 1^23 (i) are the Rabi frequencies of the 
two transitions, which may be different (because of the 
different transition dipole moments) but have the same 
time dependence since they are induced by the same laser 
pulse. 

A second example is found in rf transitions between 
the three magnetic sublevels rn = — 1, 0, 1 of a level with 
an angular momentum J = 1 in an atom trapped in a 
magnetooptical trap. The rf pulse provides the pulsed 
coupling between the m — —1 and m = sublevels, and 
also between the m = and m = 1 sublevels. The trap- 
ping magnetic field causes Zeeman shifts in the magnetic 
sublevels m = —1 and 1 in different directions but it does 
not affect the m — level [22]. This linkagepattern is 
an example of a bowtie level crossing [J, 111, la, Hi- If a 
quadratic Zeeman shift is taken into account, then the 
sublevels m = — 1 and 1 will be shifted in the same direc- 
tion, which will break the symmetry of the bowtie linkage 
and will create the "triangle" pattern of Fig. [TJ 

A third example is found in quantum rotors, for in- 
stance, in rotational ladder climbing in molecules by us- 
ing a pair of chirped ultrashort laser pulses ^]. The 
energy slope is due to the laser chirp, and the static de- 
tuning Aq arises due to the rotational energy splitting. 
If the laser pulse duration is chosen appropriately then 
only three rotational states will be coupled, with their 
energies forming the "triangle" pattern of Fig. [TJ 

The fourth example is the entanglement between two 
spin-1/2 particles interacting with two crossed magnetic 
fields, a linear field along one axis and a pulsed field along 
another axis [23|. The role of the static detuning Ao 
is played by the spin-spin coupling constant. Three of 
the four collective states form a chain, which has exactly 
the "triangle" linkage pattern of Fig. [TJ In this sys- 
tem, states tjji and -03 correspond to the product states 
II) II) and It) It), whereas state ^2 is the entangled state 

(li)IT) + IT)li))/V2. 



III. EVOLUTION MATRIX 



where R(i) is an orthogonal (because H(i) is real) trans- 
formation matrix, R^^(t) — 'R.'^{t), whose columns are 
the adiabatic states Xk {k — 1,2,3), with xi having the 
lowest energy and X3 the highest energy. As we are 
only interested in the populations at infinite times, we 
need only R(±cx)), rather than the explicit function R(t). 
R(±oo) can be easily obtained using the asymptotic be- 
havior of H(t) at infinite times. 
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(5) 



The Schrodinger equation in adiabatic basis reads 

zA(i)=HA(t)A(i), (6) 

with UA{t) = R^(t)H(t)R(t) - iR^(t)R(i), or 



UA{t) = 



Ai(t) -ivuit) -wiiit) 

-W2l{t) X2{t) -iU23{t) 
-W3l(t) -ij^32(0 h{t) 



(7) 



where the nonadiabatic coupling between the adiabatic 
states Xk{t) and x;(i) is 



m{t) ^ {xk{t)\xi{t)) ^ -'^ik{t). 



B. Assumptions 



(8) 



Our approach is based on two simplifying assumptions. 
First, we assume that appreciable transitions take place 
only between neighboring adiabatic states, Xi(t) *~* X2{t) 
and X2(0 ^^ X3(0, but not between states Xi(t) and 
X3(i), because the energies of the latter pair are split by 
the largest gap (cf. Fig. [TJ). Second, we assume that 
the nonadiabatic transitions occur instantly at the corre- 
sponding avoided crossings and the evolution is adiabatic 
elsewhere. This allows us to obtain the evolution ma- 
trix in the adiabatic basis by multiplying seven evolution 
matrices describing either LZ nonadiabatic transitions or 
adiabatic evolution. 



An exact solution of the Schrodinger equation ^ for 
the Hamiltonian ([TJ) is not known. We shall derive an 
approximation, which is most conveniently obtained in 
the adiabatic basis. 



A. Adiabatic picture 

The adiabatic states are defined as the eigenvectors 
Xk{t) (k = 1,2,3) of the instantaneous Hamiltonian 
H(i). The corresponding adiabatic amplitudes A[t) = 
[Ai{t) , A2{t) , A^it)] and the diabatic ones C{t) are re- 
lated as 



C{t) = R(t)A(t), 



(4) 



C. Evolution matrix in the adiabatic basis 

The adiabatic evolution matrix U (cxd, — cxj) is most 
conveniently determined in the adiabatic interaction rep- 
resentation, where the diagonal elements of HA(t) are 
nullified. The transformation to this basis reads 



A(t)=M(i,to)B(t), 



where 

M(t,<o) 



A 



-iAi(t,to) 

g-JA2(t,to) 

g-iA3(t:to) 

t 



k(t,to)^gXk{t')dt' 



(9) 

(,10a) 
(10b) 



and to is an arbitrary fixed time. The Schrodinger equa- 
tion in this basis reads 



iB(i) = HB(t)B(t), 



(11) 



where ^eff (0) is the effective couphng between states ipi 
and ip3 at crossing time i = 0; it is determined by the 
sphtting between the adiabatic curves A2(i) and Xsit), 



with 



HB(t) = -^ 





iA2l(t,to) 



jy^^QiJ^12{t.to) j^^^giAl3(t,to) 



1/21 e 



1^236 



iA23(*,to) 







(12) 

where Aki{t,to) = Ak{t,to) — A/(t, ip)- In this basis the 
propagator for adiabatic evolution is the identity matrix. 
The LZ transitions at the crossings at times — r, 0, r 
are described by the transition matrices [30J 



^e//(0) = A2(0) - A3(0) = i (-Ao + ^JaI + 2fi2 

(16) 
The propagator in the adiabatic basis reads 



U^(oo, -oo) = M(oo, r)ULz(T)M(T, 0)IJlz{0) 

xM(0,-T)ULz(-T)M(-r,-oo). (17) 



Ulz(-t) = 
Ulz{0) = 
Ulz{t) = 



i(p — 
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pT 
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(13a) 
(13b) 
(13c) 



D. Propagator and transition probabilities in the 
diabatic basis 



Below we present the diabatic propagator in an explicit 
form. For simplicity, we assume equal couplings 



where Pi^ (>«: = — ,0, +) is the LZ probability of nonadi- 
abatic transition and q^ is the no-transition probability 
at the crossings at times — r, 0, r, 



Here 



P>. 



q>. 



I -P.. (14) 



a_ = n,2{~T)/{2Afl\ 
ao = n,ff{0)/2A'/^, 
a+ = f}23(T)/(2A)i/2, 



(15a) 
(15b) 
(15c) 
arg r(l - tal) + J + a2^(ln a^ - 1), (15d) 



ni2{t) = n23{t) = ^{t), (18) 

although our approach is valid in the general case. This 
constraint is not applicable for the ladder climbing sys- 
tem, considered in Sec. IIIB[ where the couplings are nat- 
urally different due to the different transition dipole mo- 
ments, but is still valid for the other systems discussed. 
Then Afe;(0,-t) = Afe;(f,0), a+ = a^ = a, 0+ =(/)_. = 0, 
P+=P- = P: and q+ = q^= q. 

We find the propagator in the original diabatic ba- 
sis by using Eqs. ^, (O and pT|) as U(oo,— cx)) = 
R(cxd)U'^(oo, — cxd)R^(— cxd), or explicitly, 



U(cxD, — CX)) = 



giipi+iipa 



^-i(t>~iLpi+iip2 



pq^ 



g-2i0 _ g-i0o+2iip 



. PPO 

i(t>-iclio+itfil+iV2 /pnqo 



_g-2iipi+2jcp2p ^ g2i0-iipn+2iip2„ /gT 



Vpggo 



gi0O + 2»V3^ 

_^i4>+iip2+iV3 /-qpT 



/PPa 
(19) 

with ifi — Ai2(t, 0), ip2 = Ai2(cxD,0), ip3 — Ai3(cxD,0). The transition probability matrix, i.e. the matrix of the 
absolute squares of the elements of the propagator (I19p . reads 



PPo qpo go 

qp + pqaq + 2qpy^cos 7 p^ + q'^qo - 2qpy^cos 7 qpo 
q^ + qop^ - 2qpJq^cos 'f qp + pqoq + 2qpJq^cos 7 ppo 



(20) 



r 



where 



7 = 20- 



2^, 



(21) 



The element at the m-th row and the n-th column 
of the matrix (|20p is the transition probability Pn^m, 



that is the population of state m at infinite time, when 
the system starts in state n in the infinite past. The 
survival probabilities P1-+1 and Pa^s coincide with the 
exact expressions conjectured [3] and derived exactly for 
constant couplings [13, [3 earlier. 

In Eq. (j20|) we recognize interference terms, which arise 
because of the availability of two alternative propagating 
paths in the Hilbert space. There is also a symmetry 
with respect to the skew diagonal due to the equal cou- 
plings between neighboring states (HH]) and the equal (in 
magnitude) slopes of the energies of states i/'i and i/'s ID)- 



E. Conditions of validity 

As already stressed, our approach presumes that the 
nonadiabatic transitions occur in well-separated confined 
time intervals. This means that the characteristic transi- 
tion times are shorter than the times between the cross- 
ings, or ttransition ^ T- The transition times for diabatic 
(fi^ ^ A) and adiabatic (fi^ » A) regimes are [3l| 



''transition 
''transition 



y^2n/A, diabatic regime, (22a) 
2il/A, adiabatic regime. (22b) 



This leads to the following conditions for validity: 



Aq > VzttA, diabatic regime, (23a) 

Aq > 257, adiabatic regime. (23b) 

We shall demonstrate that the LZ-based approxima- 
tion pop outperforms its formal conditions of validity 
([23ll and is valid beyond the respective ranges. 

F. Case of Aq < and/or A<0 

Above we assumed that Aq > and A > 0. Now we 
consider the cases Aq < and A < 0. We assume that 
the couplings are even functions, n{—t) — il{t). 

Negative static detuning (Ap < 0). 

The Schrodinger equation for the propagator 
U(Ao;i,t,) is 

z|-U(Ao; i, ti) = H(Ao, i)U(Ao; t, t,). (24) 
ot 

By changing the signs of Ag , i and ti in Eq p4)) , we obtain 
the same equation, but with the f2(i) replaced by — n(t) 
[see Eq. (dJ]. It is easy to see that the change of sign of 
fl{t) is equivalent to the transformation U ^ U' = QUQ 
where Q is the diagonal matrix Q = diag{l, — 1, 1}. 
Hence we find 

i^U'(-Ao; -t, -k) = H(Ao, i)U'(-Ao; -t, -t,). 

(25) 
Because the initial condition at t — > — cx) for U(Ao;i, ii) 
and U'(— Aq; — i, —ti) at i ^ — cx) is the same, 
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FIG. 2; (Color online) The transition probabilities Pm^n 
for the transition m -^ n vs the detuning Aq for A = 
lOO/T^, flo = lO/T. Each frame compares the numerical 
(dashed red) and analytical (solid blue) results. 



we conclude that U(Ao; t, ti) — U'(— Aq; — t, —ti); hence 

U(— Ao;oo, — 00) = QU(Ao; — oo,oo)Q 

= QU(Ao;oo,-oo)^Q. (27) 

Therefore 

P™^„(-Ao) - P„^™(Ao), (m,n= 1,2,3). (28) 

Negative chirp rate {A < 0). 

We notice that Hii{A) — H33{—A), i.e. the change of 
sign of A is equivalent to exchanging the indices 1 and 
3. Hence the probabilities for A < are obtained from 
these for ^ > using the relation 

P™^„(-A) = P4_„^4-„(A), (m,n= 1,2,3). (29) 



IV. COMPARISON OF ANALYTICAL AND 
NUMERICAL RESULTS 

Below we compare our analytical approximation with 
numerical simulations. We take for definiteness the cou- 
plings in Eq. ([1]) to be Gaussians, n,{t) = fioe"* ^"^ ■ 

Figure [5] shows the nine transition probabilities vs the 
static detuning Aq. An excellent agreement is observed 
between analytics and numerics, which are barely dis- 
cernible. This agreement indicates that the dynamics 
is indeed driven by separated level-crossing transitions 
of LZ type. The analytic approximation (|20p is clearly 
valid beyond its formal range of validity, defined by con- 
ditions ((23)) . which suggest |Ao| > 25/r for the param- 
eters in this figure. The figure also demonstrates that 
the detuning can be used as a control parameter for the 
probabilities in wide ranges. 
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FIG. 3: (Color online) The transition probabilities Pm^n 
for the transition m, ^ n vs the energy slope A for Ao = 
30/r, fio = 10/r. Each frame compares the numerical 
(dashed red) and analytical (solid blue) results. 



For Ao > the five probabilities on the first row and 
the last column vary smoothly, in agreement with the 
analytic prediction. The two-photon probability Pj,^i 
vanishes rapidly with Aq, as expected, at a much faster 
pace than the other probabilities. The other four proba- 
bilities Pi^2, Pi^3, P2^2 and P2^3 exhibit oscillations, 
in agreement with the analytic prediction, due to the ex- 
istence of two alternative paths of different length from 
the initial to the final state (see Fig. [T|), with an ensu- 
ing interference. It is noteworthy that these oscillations, 
due to path interference, are not particularly pronounced, 
which might be a little surprising at first glance. How- 
ever, a more careful analysis reveals that when a control 
parameter is varied, such as the static detuning Aq here, 
it changes not only the relative phase along the two paths 
(which causes the oscillations) , but also the LZ probabil- 
ities Ptt and q^ {>c = — ,0, +)• Indeed, as Ao increases, 
we have p± — > 1 because the crossings at times ±t move 
away from the center of the pulses and 17 (±t) — > 0. These 
probabilities affect both the average value of Pm^n and 
the oscillation amplitude, with Pm^n tending eventually 
to either or 1 for large Ap, while the oscillation ampli- 
tude (which is proportional to p±) is damped. 

Similar conclusions apply to the case of Aq < because 
of the symmetry property ((28l) . It is easy to see from 
here that the survival probabilities Pn^n {n = 1,2,3) 
are symmetric vs Aq, as indeed seen in Fig. [21 

Figure [3] displays the transition probabilities vs the 
chirp rate A. An excellent agreement is again observed 
between analytics and numerics. We have verified that 
the analytic approximation ([^0]) is valid well beyond its 
formal range of validity conditions ((23)) . which suggest 
\A\ < 140/T^ for this figure; this is not shown because 
our intention here is to show the small-A range that ex- 
hibits interference patterns. As with the static detuning 



3 
O 



1.0 
0.8 
0.6 
0.4 
0.2 


1.0 
0.8 
0.6 
0.4 
0.2 


1.0 

0.8 
0.6 
0.4 
0.2 




Rabi Frequency Q. (units of 1/T) 

FIG. 4; (Color online) The transition probabilities Pm~>n 
for the transition m — *■ n vs the Rabi frequency f2o for 
Ao — 10/r, A — ZQ/T^ . Each frame compares the numerical 
(dashed red) and analytical (solid blue) results. The vertical 
dashed lines for Pi_,i, P3->3 and Pa^i show the values f2i/2 of 
the Rabi frequency for half population in the relevant states, 
predicted by our model, Eqs. pO|l and (ISlf) . 



in Fig. [21 this figure demonstrates the symmetry with 
respect to the sign inversion of A, derived in Eq. (|29p : 
the change A -^ —A \s equivalent to the exchange of 
the indices 1 and 3. The observed additional symmetry, 
P2^i = P3^2 and Pi^2 = Pj^a, is a consequence from 
the assumptions of equal Rabi frequencies and equal (in 
magnitude) slopes of the energies of states ■01 and -03 ■ 
The figure also shows that, with the exception of the 
survival probabilities P„— ,„ {n = 1,2,3), all other prob- 
abilities are asymmetric vs the chirp rate A, unlike the 
two-state level-crossing case. For A > 0, as for Aq > in 
Fig. [21 oscillations are observed in the four probabilities 
in the lower left corner but not for the probabilities in 
the top row and the right column. On the contrary, for 
A < 0, oscillations are observed only in the four prob- 
abilities in the top right corner. As discussed in regard 
to Fig. [21 the observation of these oscillations is in full 
agreement with their interpretation as resulting from in- 
terference between two different evolution paths to the 
relevant final state. 

Like the static detuning Aq, the energy slope A can be 
used as a control parameter because it affects the prob- 
abilities considerably. Around the origin {A = 0) the 
system is in adiabatic regime, while for large \A\ it is in 
diabatic regime. For instance, when the system is ini- 
tially in ^1, around the origin (A = 0) the population 
fiows mostly into state V'3, following the adiabatic state 
Xi(i). On the contrary, for large A it eventually returns 
to ipi (not visible for the chirp range in Fig. [3]). 

Diabatic and adiabatic regimes are easy to identify also 
in Fig. [31 where the nine probabilities are plotted vs the 
peak Rabi frequency Qq, which is another control param- 



eter. Consider our system initially prepared in state i/'i- 
For weak couplings the system evolves diabatically and 
therefore it is most likely to end up in the same state tpi . 
As the couplings increase, the system switches gradually 
from diabatic to adiabatic evolution; for strong couplings 
the evolution proceeds along the adiabatic state Xi(^)i 
and we observe nearly complete population transfer to 
state V'a- 

Returning to the issue of oscillations, such are barely 
seen in Fig. H) As discussed in relation to Fig. [21 a 
varying control parameter changes, besides the relative 
phase of the interfering paths, also the probabilities p^^ 
and q^, which eventually acquire their asymptotic values 
of or 1; in these limits the oscillations vanish. The 
probabilities depend on the peak Rabi frequency Qq much 
more sensitively than on the static detuning Aq and the 
energy slope A; consequently, clear oscillations are seen 
vs Aq and A, but not vs Qq, because the dependence of 
P;tc on Oq is strongest (essentially Gaussian), and hence 
the approach to the asymptotic values of the probabilities 
is fastest. 
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FIG. 5: (Color online) The final populations of states ipi, 
^2 and ^3 vs the chirp rate A for fixed Aq = lOO/T and 
f^o ~ 36.2/r, provided the system is initially in state ip^. 
The three curves cross at about A ~ 74.5/r^, indicating the 
creation of a maximally coherent superposition with popula- 
tions Pi = P2 = Pa = 1/3, which is very close to the solution 
of Eqs. (|32a)l and (|32b|) . A = 73.6/r^ shown with a vertical 
dashed line. 



V. APPLICATIONS OF ANALYTICS 

In this section we shall use our analytic approxima- 
tion for the transition probabilities ([^0]) to derive several 
useful properties of the triple-crossing system. 



A. Analytical linewidth 

We begin by deriving approximate expressions for the 
Rabi frequency required to reach 50% population in the 
n-th state for the transition m ^ n. Simple expressions 
are found for the transition 3^1, 



Oi 



'2A\n2 + AQ\/TTA\n2 



/2 



(30) 



determine by means of our analytical model values of 
Ao, A and JIq, so that we achieve arbitrary preselected 
populations at the end. For example, for a maximally 
coherent superposition state, i.e. Pi = P2 = -P3 = 3, we 
need P = 5 and po = | • This yields the following set of 
equations for Aq, rig, and A: 



ie2^oM'ln2-Ao 



7rln3/2 
A 



-2 In 3/2 = 0, (32a) 

2Aln2^At,/A\ (32b) 



An example is shown in Fig. [5] where the three final 
probabilities Pj^i, ^3^2 and ^3^3 are plotted versus 
the chirp rate A. The three probabilities cross (indicating 
the creation of a maximally coherent superposition state) 
approximately at the value predicted by Eqs. (|32ap and 
(|32bp . shown by the vertical line. 



and for the transitions 1^1 and 3 
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where a = exp (2Ao/A^). These values are indicated by 
vertical lines in Fig. |4] and are seen to be in excellent 
agreement with the exact values. 



B. Creation of superpositions 

If we prepare our system initially in state ipi and use 
A < 0, or in state V'3 and use A > 0, it is possible to 



VI. COMPARISON WITH THE EXACTLY 
SOLUBLE CARROLL-HIOE MODEL FOR Ao = 

For Ag — and constant couplings, the Hamiltonian 
^ allows for an exact solution - this is the CarroU-Hioe 
(CH) bowtie model [J|. The transition probability matrix 
for the CH model reads 



CH 



where 



pI 2^,(1 -Pc) (1-Pc)' 

2p,{l-p,) {l-2p,f 2p,(l-pe) 
{l-p,f 2p,{l-p,) pI 



(33) 

(34) 

We use this exact result as a reference for the Aq = 
limit of our approximate method, applied for constant 
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For a~^\ they read, up to 0(e '^° ), respectively 
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FIG. 6: (Color online) Comparison of the probabilities (|33[) in 
the exactly soluble Carroll-Hioe model (dashed red line) with 
our approximate solution (|20p (solid blue line) for Ao = as 
functions of the chirp rate A. Here Q, = l/T. 



e-^'^^l^ 1 - e-'^"'/^ 

1 



(36a) 



CH 



2e"™'/2 1 - 2e~''°'/2 

1 - 2e-'^"'/2 2e-'^"'/2 



Equations (|35ap and (|35bp demonstrate that our ap- 
proximate solution (j20p reproduces well, for some prob- 
abilities even exactly, the correct small-a asymptotics, 
which corresponds to the large-A ranges in Fig. [Sj The 
reason is that the small-a (diabatic) regime corresponds 
to weak coupling; in the pcrturbative regime the pres- 
ence of level crossings, let alone their distribution in time, 
is less significant, fn the large-a (adiabatic) regime the 
crossings become very important and definitive for the 
dynamics. Then Eq. (|36ap deviates from the correct 
asymptotics (j36bp . but still has the correct asymptotic 
values for a — > cxd. The correct, or nearly correct, small- 
a and large-a asymptotics of our approximate solution 
(j20p explain its surprising overall accuracy in Fig. [6l 



coupling r2(i) = 51 = const. We emphasize that tak- 
ing this limit is an abuse of the method because in the 
derivation we have assumed that the crossings are sepa- 
rated, which has justified the multiplication of propaga- 
tors. Nonetheless, it is curious and instructive to push 
our approximation to this limit. For Aq = the LZ pa- 
rameters are a± — n/\/2A — a and ao = a/2. Therefore 
we have Pq ~ p^ = p. 

Figure [6] presents a comparison between the exact 
Carroll-Hioe solution p3p and our approximate solution 
(|20p . Quite astonishingly, our approximate solution is 
not only qualitatively correct but it is even in a very 
good quantitative agreement with the exact solution; 
we witness here yet another LZ surprise where our LZ- 
based model outperforms expectations in a limit where 
it should not be adequate. 

The observed feature of our approximate solution can 
be explained by examining the asymptotics of the approx- 
imate probabilities (PD|) and the exact CH values ([55)1 for 
a <C 1 and a 3> 1 . For a ^ 1 the approximation ((20)) and 
the CH solution ([55)1 read, up to 0(a^), respectively 



CH 



— 57ra^/4 7ra^ 7ra^/4 






na^ 1 — 27ra^ na^ 


,(35a) 


7raV4 7ra2 1 - bira^/A 




1 — na^ -KO^ 






■no? 1 — 27ra^ ■no? 




(35b) 


7ra 1 — 7ra 







VII. DISCUSSION AND CONCLUSIONS 

We have developed an approximate analytical model 
that describes the time-dependent dynamics of a quan- 
tum system with three states, which have linearly chang- 
ing energies of different slopes and are coupled with pulse- 
shaped interactions. Our approach is based upon the 
two-state LZ model, i.e. we assume independent pairwise 
transitions between neighboring states, described by the 
LZ model. We have performed detailed comparison of 
our analytic approximation with numerical simulations, 
versus all possible interaction parameters and for all nine 
transition probabilities, which has revealed a remarkable 
accuracy, not only in smooth features, but also in describ- 
ing detailed interference features. This accuracy shows 
that indeed, the physical mechanism of the three-state 
dynamics is dominated by separated pairwise LZ transi- 
tions, even when the crossings are too close to each other. 

We have derived the formal conditions of validity of 
our LZ approach, Eqs. (P5|) . using the concept of tran- 
sition time. However, a comparison with numeric simu- 
lations has revealed that our approximation is valid well 
beyond the formal ranges of validity. One of the reasons 
is that for two of the survival probabilities, Pi^i and 
P3^3, our LZ approximation produces the exact results. 
We have found that even in the extreme case of vanishing 
static detuning, where our approach should not be valid 
because the three crossings coalesce into a triply degen- 
erate bowtie single crossing, it still produces remarkably 
accurate results because of nearly correct asymptotic be- 
haviors of the transition probabilities. 

One of the useful and interesting features of the "tri- 
angle" linkage pattern (Fig. [1]) is the presence of intrinsic 



interference effects. Our "sandwich" approach, with its 
implementation in the adiabatic interaction representa- 
tion, allows for an easy incorporation of different evo- 
lution paths in Hilbert space between a particular pair 
of states. Such path interferences are identified in only 
four of the nine probabilities. Another source of inter- 
ferences could be nonadiabatic transitions in the wings 
of the Gaussian pulses, where the nonadiabatic couplings 
possess local maxima; these interferences would be vis- 
ible in all nine probabilities. We have found, however, 
that only the path interferences are clearly identified. 

A substantial contribution to the path interferences 
is played by the LZ phases (j)^. The LZ phase is of- 
ten neglected in applications of the LZ model to multi- 
ple crossings, in the so-called "independent crossing" ap- 
proximation, where only probabilities are accounted for. 
Although such an approach occasionally works, miracu- 
lously, as in the exactly soluble Demkov-Osherov [13] and 
Demkov-Ostrovsky j^] models, the present simple, but 
very instructive model, demonstrates that in general, the 
LZ phase, as well the dynamical adiabatic phases, has to 
be properly accounted for, which is achieved best in an 
evolution-matrix approach, preferably in the adiabatic- 
interaction representation |32l | . 

In order to be closer to experimental reality, in the ex- 
amples we have assumed pulsed interactions, specifically 
of Gaussian time dependence. This proved to be no hin- 
drance for the accuracy of the model, which is remarkable 
because we have applied the LZ model (which presumes 
constant couplings) at crossings (the first and the last 
ones) situated at the wings of the Gaussian-shaped cou- 
plings where the latter change rapidly. This robustness 



of the approach can be traced to the use of the adiabatic 
basis where the pulse-shape details are accounted for in 
the adiabatic phases. 

We have used the analytic results to derive some use- 
ful features of the dynamics, for instance, we have found 
explicitly the parameter values for which certain prob- 
abilities reach the 50% level, and for which a maxi- 
mally coherent superposition is created of all three states 
Pl=P2 = P3^ 1/3. 

In the specific derivations we have assumed for the sake 
of simplicity equal couplings for the two transitions and 
slopes of different signs but equal magnitudes for two of 
the energies. These assumptions simplify considerably 
the ensuing expressions for the probabilities; moreover, 
they are actually present in some important applications 
(quantum rotors, Zeeman sublevels in magnetic field and 
spin-spin entanglement). The formalism is readily ex- 
tended to the general case, of unequal couplings and dif- 
ferent slopes, and we have verified that the resulting LZ- 
based approximation is very accurate again. 

To conclude, the present work demonstrates that, once 
again, the LZ model outperforms expectations when ap- 
plied to multistate dynamics, with multiple level cross- 
ings and a multitude of evolution paths. 
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